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ABSTRACT 

We study possible correlations between ultrahigh energy cosmic rays (UHECRs), observed by Auger, 
AGASA, and Yakutsk, and nearby active galactic nuclei (AGNs) and Fermi sources. We consider the de- 
flection effects by a Galactic magnetic field (GMF) model constrained by the most updated measurements. 
We found that the average deflection angles of UHECRs by the Galactic magnetic fields are less than 4°. A 
correlation between the Auger cosmic -ray events and nearby AGNs with a significance level of ^ 4cr was 
found for the Auger UHECR data sets with or without deflection correction. No correlation was found be- 
tween the AGASA/Yakutsk events with nearby AGNs. Marginal correlations between the Auger events and 
the Fermi sources, and between AGASA events and Fermi AGNs were found when the deflections cal- 
culated by the GMF model were considered. However, no correlation was found between the Yakutsk data 
and Fermi sources. Some Fermi sources are close to the arrival directions of UHECR events detected by 
Auger, AGASA, and Yakutsk, most of which are probably chance coincidence rather than objects producing 
UHECRs in the nearby Universe. Four Fermi sources, NGC 4945, ESO 323-G77, NGC 6951, and Cen A, 
within 100 Mpc have UHECR events within 3°1 from their positions, which could potentially be cosmic -ray 
accelerators. However, the association can only be confirmed if more UHECRs are preferably detected in these 
directions 

Subject headings: (ISM:) cosmic rays - galaxies: active - magnetic field - methods: statistical 



1. INTRODUCTION 

The spectrum, origin, and composition of ultrahigh en- 
ergy cosmic rays (UHECRs) with energies ^ 10^^ eV (=10 
EeV) a r e a lon g s tanding m ystery in high-energy astrophysics 
(iHiUas! 11984- iGreisenI ([1966) and Zatsepin & K uz'min 
(11966) showed a theoretical distant limit for the cosmic rays 
with energies of order lO'^^ eV traveling through the mi- 
crowave background radiation field, which is called the GZK 
effect. Because of the GZK effect, particles with energies 
above 10 EeV are able to reach our Earth only from nearby 
sources within about 100 Mpc. Another barrier in the in- 
vestigation of the UHECR origin is the deflections of UHE- 
CRs by the magnetic fields. Due to the poor knowledge of 
the extragalactic and intergalactic magnetic fie lds, the de- 
flections of UHECRs have not yet understood. iDolag etal] 
(|2004, 2005) suggested that the deflections by ex tragalactic 
magnetic fi elds are generally less than 1°, while iRyu et alJ 
dlOlO) and ISid et alj ( 120031) claimed that could be larger 
than 10°. The Galactic magnetic fields (GMFs) are rela- 
tively better known (e.g.lian et al. 2006; Sun et al. 200_8) and 
are widely discussed in the studies of UHECR origin (e.g. 
StaneV '1997^, 'Tinvakov & TkacheV '2002"; 'Prouza & Smidal 
S)03; Nagar & Matulich 2009). KachelrieB et al. (2007|) con- 
cluded that the deflections of UHECRs by the GMFs can- 
not be neglected even for the protons of E = lO^^eV, since 
the deflection angles are com parable with the angular r esolu- 
tion of current experiments. Nagar & MatulichI (l2009i) tried 
seven GMF models to study the correlations between UHE- 
CRs and source population(s). However, no halo compo- 
nent was included in the four GMF mod els they used and 
another three GMF models adopted from ISun et al.1 (|2008|) 
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have a strong halo component about 7 fiG. Observational con- 
straints on the Galactic magnetic field streng th ( Han & Oiaol 
ll994tlHan et"anil999HMorris & Serabvn|[T99 6) and th e con- 
figuration of disk magnetic fields dHan et al E006;.Hani2009h 
should be carefully considered in the GMF model. 

Since the discovery of UHECRs (iLinslevI Il963h . many 
equipments have b een used to sear ch for these events, in- 
cluding Fly's Ev e (iBird et al.l|1994l). Yakutsk Extensive Ai r 
Showers Array dlvanov et al. I I2003al; l^avdin et al.' "2005^ 
Akeno Giant Air Shower A rray (AGASA; [Havashida et al] 
'2000^, iTakeda et al.1 l2003h . High Resolution Fly's Eye 
cosmic-ray detector (HiRes; lAbbasietaDl2004 .20 08) and 
Pierre Auger Observatory (PAO; [Abraham et all 120041: 
iThe Pierre Auger CoUaborationll2007l l2008h . The"existence 
of the GZK cutoff has been observed by the HiRes and 
Auge r dAbbasi et al.1 l2008t iThe Pierre Auger Collaboration! 
I2OO8I) . Some objects have been suggest ed to be 
possib le sources of UHECRs, e.g., pulsars (iBlasi et al.l 

I2OOOI) . active galactic nuclei (AGN s) and subclasses 

of A GNs (Protheroe & Szabol 119921; iFar rar & Biermannl 
1298'; l^yakov & Tkachev '2001 a"b'; 'Virmani et al. 2003 
Gorbunov et al. 2002, 2004; Abbasi et al. 2006; Farrar et al 
20091) . radio lobes of FR II galaxies dRachen & Biermann 



20091), 
1991 



T99I ^Hardcast le et al.l 
1995; Milgrom&Uso3 



2009), 



1995) 



bursts IWa 



and 7-ray 

However, the real sources of 
UHECRs are not known yet. AGNs are favored as the most 
probable sources for accelerating particles to the extreme en- 
ergies (IHillasl[T984 for a long time. 

Recently, iThe Pierre Auger Coflaborationl (l2007l l2008l) 
studied the correlation between the arrival directions of 
UHECRs and the position s of nearby AGNs in the 
IVeron-Cetty & VeronI (l2006h AGN catalog (hereafter VCV 
catalog). They concluded that the arrival directions of cos- 
mic rays with energies above ^ 60 EeV are anisotropic 
and UHECRs have a good correlation with the positions of 
nearby AGNs (z < 0.018). The intriguing result attracted 
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much attention. I lvanovl ( |2008|) found the correlation between 
Yakutsk UH ECRs and the nearby VCV AGNs (< 100 Mpc). 
iGeorge eTal. (2008) investigated the correlation between the 
Swift Burst Alert Telescope AGN catalog with the Auger 
UHECR events, and found a correlation at a significance level 
of 98% when the AGNs were weighted by their hard X-ray 
flux and the Auger experiment exposure. However, some as- 
sociated AGNs of Auger events may not h ave enough energy 
to acc elerate particles to ultrahigh energies jMoskalenko et al.l 
|2009|) . The High Resolution Flys Eye Collaboration searched 
for possible correlation between the HiRes UHECRs and 
AGNs located in the northern hemisphere; however, no sig- 
nificant correlation was found. 

The Fermi high energy 7-ray sources are also possible 
UHECR sources. The recently released Fermi Large Area 
Telescope First Source Catalo g (IFGL) contains 1451 7-ray 
point sources ( Abdo et al. 2010) with nearly uniform sky cov- 
erage (Atwood et al. 2009). Mirabal & Oya (2010) first in- 
vestigated the correlation between Auger UHECRs and IFGL 
sources without considering the deflection by the GMFs and 
redshifts of Fermi AGNs, and concluded that the UHECRs 
are not associated with Fermi sources. The possible correla- 
tion of UHECRs and Fermi sources should be re-examined 
after the GZK cutoff and the UHECR deflection by the GMFs 
are considered. 

In this work, we construct a new GMF model based on the 
updated measurements of the Galactic magnetic fields and in- 
vestigate the deflections of UHECRs by the GMFs. Consider- 
ing the GZK cutoff and the deflection correction through our 
GMF model, we re-examine the possible correlation between 
UHECRs and nearby AGNs and Fermi sources. In Section 

2, we discuss available data of UHECRs detected by Auger, 
AGASA, and Yakutsk and possible astrophysical objects. The 
deflections of UHECRs by the GMFs are discussed in Section 

3. The correlation studies are given in Section 4. Discussions 
and conclusions are presented in Section 5. 

2. UHECR DATA AND POTENTIAL COSMIC RAY SOURCES 
2.1. UHECR events 

High-quality UHECR data and reasonable deflection cor- 
rection by the GMF models are crucial to understand the ori- 
gin of UHECRs. Therefore, we only consider UHECR events 
which satisfy two criteria: (1) good angular and energy res- 
olutions; (2) ultrahigh energy {E > 40 EeV) which has a 
predictable small deflection angle. 

We work on the UHECR events recorded by Auger, 
AGASA, and Yakutsk. The HiRes UHECRs have a typ- 
ical angular resolution of 0°6 (Abbasi et al. 2004) and 27 
events have been published. However, we did not use 
them since the detail position s and energies are not available 
(iKachelrieB & Semikozll20()6l and references therein). 

Auger is located in Argentina and began to collect data from 
2004 January 1 . It has two systems, one to measure fluores- 
cence in the atmosphere and the other to detect Cerenkov light 
from relati vistic partic les. The angular resolution of Auger is 
about 0°9 (lAvdlioOTh . Eighty-one events with > 40 EeV 
have been recorded by Auger from 2004 January 1 to 2007 
August 31 but only 27 events with energies above 57 EeV 
were published. Recentl y, another 31 new events (E > 5 7 
EeV) were detected (The Pierre Auger Collaboration! |20()9|) ; 
however, data are not yet available. AGASA and Yakutsk are 
all located in the northern hemisphere. AGASA has been op- 
erated for 12 years, and ceased oper ation on 2004 Januar y 
4 Its angular resolution is about 1°8 (iHavashida et al.ll2000l) . 



Up to 2000, 57 events with E > 40 EeV have been pub- 
lished ("Hayash ida et alJ l2000'). Yakutsk collected 51 events 
with E > 40 EeV, its angular resolution is smaller than 5° 
dPravdin et alj|2005h . Therefore, the UHECRs data used in 
this work include 135 events, 57 recorded by AGASA and 51 
by Yakutsk in the northern hemisphere with E > AO EeV, and 
27 events recorded by Auger in the southern hemisphere with 
^ 57 EeV. 

Due to the different angular resolution, different energy cal- 
ibration, and different sky exposure for Auger, AGASA, and 
Yakutsk UHECR events, we search for the possible correla- 
tions separately between the three sets of UHECRs and astro- 
physical objects. 

2.2. Potential Astrophysical Objects as Cosmic-Ray Sources 

The nearby AGNs use d in this work are extracted f rom the 
13th VCV AGN catalog ( Veron-Cettv & VeronllloTob . with a 
redshift limit of z ^ 0.024 which corresponds to a GZK cutoff 

100 Mpc for the Hubble constant i/o = 71kms~^ Mpc"^ 
The VCV catalog includes all known AGNs reported in the 
literature. There are 133,336 quasars, 1374 BL Lac objects, 
and 34,231 active galaxies. This catalog is not complete and 
not uniform due to different selection criteria and telescope 
time devoted to different sky areas. On the other hand, the 
difficulties in the classification of galaxies with weak AGN- 
like activity result in the confusion in identification of such 
sources. Thus, there are some flaws when using some given 
particular AGN cat alog to study their possi ble correlations 
with UHECRs (see iMoskalenko etani2009l) . AGNs are fa- 
vored as the most probable sources for accelerating particles 
to the extreme energies (Hillas 1984) and almost all known 
AGNs are listed in the catalog. In this work, we work on 
the VCV nearby AGNs and compare t he results with that by 
iThe Pierre Auger Collaboration! ( l2007h . For the VCV AGN 
catalog, the incompleteness is particularly serious around the 
galactic plane of |6| 10°. Therefore, only the 830 nearby 
AGNs with \b\ > 10° are used in this work. 

Another type of UHECR s ource candidates investigated in 
this w ork are Fermi sources. iThe Fermi/LAT CoUaborationJ 
(l2009l) presented the initial 3 month results for the 205 
most significant 7-ray sources with energies above 100 MeV 
Recently, after one-year observations, Fermi/LAT r eleased 
the c atalog for the all-sky 1451 7-ray sources (Ab do et alJ 
I2OIOI) . which contains 820 identified sources and 63 1 uniden- 
tified sources. The Fermi/LAT first-year 7-ray source cat- 
alog (IFGL) is a complete 7-ray source sample with a 7- 
ray flux threshold of > 4xl0~^° cm~^ s^^ in the energy 
range 1-100 Ge V.I Abdo et al.l (12010) identified many Fermi 
AGNs with redshift data. Due to the GZK cutoff, we restrict 
the Fermi sources with z ^ 0.024. There are only eight 
AGNs that satisfy this criterion. In this work, we adopt these 
eight objects and other Fermi sources without redshift data 
to search for their possible correlation with UHECRs. Be- 
cause the diffuse emission dominates at low Galactic latitude 
(|6| ^ 10°), the Fermi sources of |6| ^ 10° are also dis- 
carded. In total, 635 Fermi sources are used in this work, 
including 262 identified sources (235 AGNs) and 373 uniden- 
tified sources. 

3. THE GALACTIC MAGNETIC FIELDS AND DEFLECTION OF 
UHECRS 

The GMFs have large-scale regular and small-scale turbu- 
lent components. The deflection angles of UHECRs caused 
by the turbulent fields are typically 1 order of magnitude 
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small er than that by the regular fields dTinvakov & TkachevI 
I2005h : therefore, we ignore the turbulent component in this 
work. The Galactic magnetic fields in general are described 
as the regular magnetic field in the disk and the possible 
large-scale field in the halo. Previous magnetic field mod- 
els for the disk are either axissymmetr ic (e.g. Stanev 1997) 
or bi-symmetric (e.g. lHarari et al.l ll999; Tinv akov & TkachevI 
12002) . However, no ne of these simp l e models agree s well 
with the observations dMen et al.ll2008HSun et al.ll2008h . 
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Fig. 1 . — Configuration of the disk magnetic field. Tlie shaded area outhnes 
spiral arms given by Hou et al. ( 2009). 



Here we develope d a toy model bas ed on rotation measures 
(RMs) of pulsars bvl Han et all (l2006h . Magnetic fields in the 
disk are reversed from arms to inter-arms which has recently 
been verified by Nota & Katgert (2010). The radial profile of 
the field strength can be described as B{R) — Bq exp(— (i? — 
Ro)/Rb), where Bq — 2.1 /iG is the local field strength, 
R is the distance from the Galactic center, Rq — 8.5 kpc is 
the galactocentric distance of the Sun, and the scale radius is 
Rb — 8.5 kpc. Here, we use the four-arm model of Hou et al] 
(120091) to describe the spiral structure of our Milky Way. In 
the polar coordinates, the ith arm can be described as i? = 
Ri ex.p[{9 — 9i) ta,n.tjji], where Ri is the initial radius, 9i is the 
start azimuthal angle, and 4'i is the pitch angle of the arm. The 
va lues of t hese parameters of four arms are given in Table 1 
of iHou etal](l2()09i) . To ease the model description, magnetic 
fields within 4.6 kpc do not have reversals. The initial width 
of each arm is set to be 0.4 kpc in our work. The pitch an gle of 
the magnetic field is —11° as used bv lHan et alj (l200d) . The 
configuration of the disk magnetic fields in our toy model is 
displayed in Figure [T] which gives the counterclockwise field 
in the arms and the clockwise field in the inter-arm regions. 
Here, we remind that this is only a toy model and the four- 
arm model is not the best one to match all tracers for the spiral 
arms (Hou et al. 2009). 

The halo magnetic fields consist of a dipole poloidal field 
and a toroidal field with opposite directions above and be- 
low the Galactic plane. The field configuration was derived 
from the antisym metric RM sky revealed by the extragalac- 
tic radio sources dHan et alJ [T997. 199 9) and the vertical fila - 
ments in the Galactic center lYusef-Zadehet"anil984 12004 . 



At present, it is not clear whether such toroidal fields extend 
from the solar vicinity to the Galactic center and what the 
scaleheight they have. Here we use the formula given by 
iProuza & Smidal d2003h in the Cartesian coordinates to de- 
scribe the toy model: 



Bx = -Bt sin((/))sign(z) 
By = Bt cos((^)sign(z), 



(1) 



where Bt is given by 



Bt — 



Bn 



^/aF+y^ < R 



1 



■ exp(- 



R 



-), y^x^ +y^> R 
(2) 

Here i? = 15 kpc is the radius of the toroidal field, H — 
1.5 kpc is the scaleheight of the toroidal disk, P = 0.3 kpc is 
the halfwidth of the Lorentzian distribution, and the maximal 
value of the magnetic field is iJmax = 1 A*G- 

The poloidal field is modeled as a dipole with a cylin- 
der (height 300 pc, diameter 200 pc) in the Gal actic center. 
The ex isten ce of the dipole fie ld is questioned bv lTavlor et al.l 
d2009l) and iMao et al.l (l2010h . However, their data in the 
Galactic pole region clearly support the local vertical field. In 
this work, we try to include the dipole field in our toy model. 
In the polar coordinates the poloidal field strength is (in the 
zz— plane, 9 ranges from to tt and goes from north to south 
pole) 



B = Bp \/3 0082(6*) + 1. (3) 
Then the Cartesian components of the poloidal field are 



where 



Bp = 



lOO/Rl 
1 

0.2/R% 
2000 



l.5Bp sin 20 cos 
1.5Bp sin 29 sin cj) 
Bp{3 cos^ 9-1), 



5 < i?p < 15 

2 < i?p < 5 

0.1 <Rp <2 

Rp < 0.1 and \z\ < 0.15. 



(4) 



(5) 



Here Rp 



z^. Rp and z are in units of kpc, and Bp 
is in units of /iG. The constants in Bp were selected to meet 
the characters of observed filaments in the Galactic center (1- 
2 mG; Morris & Serabyn 1996), and a 0.2 /iG vertical field 
component in the vicinity of the Sun dHan & Qiao 1119941) . 

In our toy model, the field transition from the arms to the 
inter-arms is not smooth, the influence of the bar in the Galac- 
tic center is not considered yet. The detailed model is beyond 
the scope of this paper We will investigate how much UHE- 
CRs are reflected by each magnetic field component in the 
model. 

The cosmic rays are d eflected in the GMFs be cause of the 
Lorentz force. FoUowing' KachelrieB et alj d2007h . the net de- 
flection can be approximated as 







Zqe 
pc 



dl X Bf 



(6) 



where Zq^ is the charge of cosmic -ray particles, p is the 
momentum along the line of sight (LOS), and Bt is the field 
component perpendicular to the LOS. The integral is along 
the LOS from the source to the observer. In our work, the 



Fig. 2. — Deflection maps of CR protons of various GMF components and their combinations in our toy model and tlie PS model jProuza & Smfdal2003l) . The 
energy of CRs is fixed to be 40 EeV. All the maps are plotted in the Galactic coordinates. Each arrow goes from the observed direction of a CR on the Earth to 
the actual incoming direction of the source. 



Hammurabi code dWaelkens et alj|2b09[) is used to calculate 
the all-sky deflections. 

In Figure|2l we show the deflection angle maps of UHECRs 
for various magnetic field components and their combinations 
in the toy model, and also the PS model ( Prouza & Smidal 
l2003h which has a different disk component. The energy of 
CR protons is fixed to be 40 EeV. The disk field generates a 
strong deflection near the Galactic disk and in the high lati- 
tude. The toroidal field is similar to the disk field and has a 
large deflection near the disk, while the orientations are oppo- 
site above and below the disk. The dipole component has very 
strong deflections in the Galactic center region. Deflection 
maps by our toy model and the PS model have some similar 



features: (1) the maximal deflections take place in the inner 
Galactic disk and Galactic central regions; (2) in the south- 
ern sky of the outer Galaxy, the deflection is as strong as that 
in the Galactic center, while in the northern sky, a very small 
deflection occurs. This is a result of combining the toroidal 
and the dipole magnetic field components. However, there is 
quite difference in some regions between these two models: 
in the northern sky, our model generates smaller deflections 
from ^ = 0° to ^ = 180°, and the orientation is contrary to 
the PS model in the region of -90° < / < -180°; in the 
southern sky, the PS model has larger deflections, particularly 
in the region of -90° < I < -180°. 

The energy of observed UHECRs in our data sets varies 
from 40 EeV to more than 200 EeV. The deflections for the de- 



Fig. 3. — Deflection of 135 detected UHECRs according to our GMF toy model (in the Galactic coordinates). The circles ai'e the original airival directions 
of 135 UHECRs observed on the Earth, the squares, triangles, and dots are the deflection corrected positions of the Auger, AGASA, and Yakutsk CR events, 
respectively. The size of symbols is proportional to the CR energy. The star marks the Cen A. 
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Fig. 4. — Average deflection angles of 135 detected UHECRs as a function 
of energy by our GMF toy model and the PS model (Prouza & Smfda 2003). 
The solid lines are our model (black dot) and PS model (black diamond), 
dotted lines are three GMF components, dashed lines are their combinations. 



tected UHECR events according to our toy model are shown 
in Figure |3] Figure |4] shows the average deflection angles of 
135 UHECRs by the two GMF models and the various mag- 
netic component combinations. We found that the deflections 
decrease with the CR energy as expected. Magnetic fields in 
our model generate smaller deflections than those in the PS 
model. CR events with energy below 90 EeV are deflected 
by an angle < 3°. The arrival directions of the events with 
energy above 90 EeV are deflected less than 2°. The deflec- 
tion angles are generally similar to the angular resolution of 
the UHECR detectors ( - 1° - 2°); therefore the GMF de- 
flection correction is still desired in understanding the origin 
of UHECRs. 

4. CORRELATIONS BETWEEN THE UHECR EVENTS AND 
POTENTIAL COSMIC-RAY SOURCES 



In this work, we try to correct the deflections of the UHECR 
events by the GMFs, and then search for possible correla- 
tions between three UHECR data sets (Auger, AGASA, and 
Yakutsk) and nearby AGNs and Fermi sources. We will com- 
pare the correlation results with and without the deflection 
correction. Here, we first introduce the correlation analysis 
method. 



4. 1. The Correlation Analysis Method 

To carry out the correlation analysis, we use the 

angular correlation functio n method descri bed in 

Tinvakov & TkachevI ^2W1^l iGorbunov etaH (l2002i 
2004 . and iGorbunov & TroitskyI (l2005h . For a sample 
with Ur cosmic -ray events, we count the number Ncs for 



UHECR-source pairs within a given angle 6, which is called 
"bin size" and various from to a large angle. We count 1 
if at least one potential cosmic-ray source (such as AGNs) 
falls into "the bin" and count if no source falls into the 
bin. To check the chance probability, we use Monte Carlo 
simulations. We generate a large number (e.g., with N= lO"', 
10^ for some cases) of simulated sets of UHECR events, 
each set has the same number of events as the real sample. 
The simulated UHECRs are isotropic and the locations of 
simulated events are random. The distribution of locations is 
constrained by the overall exposure of the UHECR detectors. 
For a given S, we first count the number of UHECR-source 
pairs N„ic for each simulated UHECR sample, and obtain 
a mean Nmc and the variance (j,nc from N simulation sets. 
The exposure of Auger, A GASA, and Yakutsk depends on 
the c elestial declination jThe Pierre Au^er Collaboration! 
l2008t ImedaetaH [19991: llvanov et all I2003bh . The Auger 
exposure used in our work is obtained from the fitting of 
the declination distribution of low-energy Auger events 
(E < 10 EeV, from the Auger Web site^) with a third-order 
poly nomial function. The AGASA exposure function is taken 
from iTakedaetalJ (I1999I). T he Yakutsk exposure function 
is derived from llvanov et al.l (f2003h) . The probability that 
the observed UHECR-source pairs are in coincidence with a 
random distribution is estimated by 

' http://auger.colostate.edu/ED/ 



6 



TABLE 1 

Correlation Between UHECRs with AGNs and Fermi 7-Ray Sources With or Without Deflection Corrections by GMFs 







Auger Events 


AGASA Events 


Yakutsk Events 


(1) 


(2) 


(3) (4) (5) (6) 


(7) (8) (9) (10) (11) 


(12) (13) (14) (15) (16) (17) 


Sources 


GMF 


SONcsNrrrc P 


a <S(°)iVcsiVmc P 


a SONcsNrnc P a 



VCV/AGNNon 2.4 17 

Toy model 3.6 20 

Disk 3.3 21 

Dipole 2.1 14 

Toroidal 3.3 20 

Disk+Dip. 3.0 18 

Disk+Toro. 3.9 21 

Toro.+Dip. 3.3 19 

PS model 3.0 18 

IFGL/All Non 4.0 17 

Toy model 0.9 5 

Disk 2.4 10 

Dipole 4.2 19 

Toroidal 0.6 2 

Disk+Dip. 1.8 9 

Disk+Toro. 1.8 7 

Toro.+Dip. 3.6 17 

PS model 2.1 12 

IFGL/AGNNon 7.6 20 

Toy model 0.9 2 

Disk 8.1 20 

Dipole 3.3 9 

Toroidal 7.5 19 

Disk+Dip. 1.8 4 

Disk+Toro. 8.1 20 

Toro.+Dip. 3.9 9 

PS model 5.4 13 

IFGL/Un-idNon 7.0 21 

Toy model 3.0 11 

Disk 6.0 19 

Dipole 1.2 4 

Toroidal 8.1 21 

Disk+Dip. 0.9 3 

Disk+Toro. 7.8 21 

Toro.+Dip. 3.6 12 

PS model 1.8 8 
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2x10"^ 
2x10"^ 

0.091 
2.5x10"= 
0.056 
0.024 
0.073 
5.3x10"^ 



4.6 
4.2 
4.2 
4.2 
4.2 
1.5 
4.1 
1.9 
2.1 
2.4 
3.0 



4.0 9.3x10" 



13.8 

1.2 

15.1 

5.1 

8.7 

18.3 

6.1 

15.3 

1.2 

19.0 

0.7 

18.6 

8.1 

2.5 



0.012 
0.030 
0.010 
0.045 
0.047 
0.100 
0.021 
0.051 
0.031 
0.120 
0.030 
0.085 
0.066 
1.7x10" 



•^2.8 
2.4 
2.7 
2.4 
2.0 
1.9 
1.5 
2.4 
1.8 
2.6 
1.5 
2.9 
1.6 
1.7 

^ 3.7 



1.4 
2.1 
0.6 
2.7 
3.6 
1.5 
2.4 
2.4 
1.5 
3.0 
2.4 
2.7 
5.7 
4.8 
1.8 
3.6 
5.4 
1.5 
3.8 
2.7 
2.7 
4.5 
4.2 
1.8 
3.6 
1.8 
1.5 
6.0 
5.4 
5.4 
6.0 
6.3 
6.3 
6.0 
5.1 
0.6 



11 6.7 

17 11.5 
4 1.4 
22 18.3 
29 25.1 

13 7.1 

18 14.6 
15 14.6 
11 6.7 
25 20.9 
18 12.9 

22 17.1 
45 42.0 

41 35.2 

14 8.4 
31 24.9 

42 38.0 
11 5.6 

23 15.8 

15 7.8 7, 
15 8.6 
27 20.0 

17.4 
4.0 5, 
13.4 
3.9 
2.6 7, 



20 
10 
20 
6 
9 



36 32.6 
29 25.6 
32 27.8 

37 31.9 
35 32.2 

38 32.2 
32 30.0 
29 24.9 

1 0.5 



0.066 

0.049 

0.048 

0.172 

0.164 

0.019 

0.179 

0.499 

0.062 

0.153 

0.070 

0.095 

0.180 

0.038 

0.033 

0.056 

0.106 

0.019 

0.025 

0x10" 

0.020 

0.032 

0.255 

2x10" 

0.027 

0.183 

0x10" 

0.201 

0.196 

0.151 

0.093 

0.251 

0.052 

0.335 

0.151 

0.387 



1.8 
1.9 
2.3 
1.1 
1.1 
2.4 
1.1 
0.1 
1.8 
1.2 
1.7 
1.5 
1.1 
1.9 
2.1 
1.7 
1.4 
2.4 
2.2 

3 2.8 
2.4 
2.0 
0.8 

^ 3.2 
2.1 
1.1 

U.l 
1.0 
1.0 
1.2 
1.5 
0.8 
1.7 
0.6 
1.2 
0.7 



3.0 
3.6 
0.3 
3.3 
3.0 
2.4 
3.0 
3.6 
1.2 
2.6 
4.2 
4.2 
4.2 
5.1 
3.6 
5.1 
3.9 
3.3 
5.2 
3.6 
0.6 
4.2 
3.6 
3.6 
3.9 
3.9 
3.3 
2.6 
6.6 
5.1 
6.9 
5.1 
6.0 
5.7 
7.5 
6.6 



19 
25 
1 

22 
20 
15 
21 
25 

7 
16 
28 
29 
29 
32 
25 
33 
27 
23 
21 
15 

1 

17 
12 
16 
15 
15 
12 
10 
30 
23 
30 
23 
28 
28 
32 
30 



17.9 
21.4 
0.3 
19.9 
17.4 
13.4 
17.4 
21.4 
4.2 
11.7 
22.7 
23.9 
23.3 
27.7 
19.7 
27.7 
20.7 
16.9 
19.8 
11.6 
0.4 
14.9 
11.6 
12.2 
13.1 
13.1 
10.3 
5.5 
23.3 
17.4 
25.2 
16.6 
21.6 
19.3 
26.8 
23.9 



0.420 
0.155 
0.277 
0.307 
0.237 
0.350 
0.148 
0.155 
0.117 
0.096 
0.053 
0.063 
0.042 
0.080 
0.066 
0.034 
0.028 
0.042 
0.416 
0.150 
0.343 
0.293 
0.498 
0.124 
0.305 
0.305 
0.312 
0.038 
0.015 
0.052 
0.072 
0.029 
0.029 
2.9x10" 
0.041 
0.028 



0.4 
1.2 
1.2 
0.7 
0.9 
0.5 
1.2 
1.2 
1.5 
1.5 
1.8 
1.7 
1.9 
1.6 
1.7 
2.0 
2.1 
1.9 
0.4 
1.2 
0.9 
0.7 
0.2 
1.3 
0.7 
0.7 
0.6 
2.1 
2.2 
1.8 
1.6 
2.1 
2.0 
2.8 
1.9 
2.0 



number of simulated sets with N„ic ^ Nc, 
N 



(7) 



The significance of the correlation can be defined as a{6) = 
N^s-N^c _ emphasize that P and <7 vary with 5. Larger 
chance probability P{5) indicates that the observed pairs of 
cosmic rays and astrophysical objects are more likely the 
statistical coincidence of random isotropic UHECR events. 
The higher significance corresponds to a smaller P{5) which 
suggests that the objects of the pairs are more likely to be 
potential UHECR sources. In this paper, P < 10"^ is 
believed as an indicator of some correlation. We use this 
method to evaluate the possible correlations between the 
Auger/ AGASA/Yakutsk events and the possible cosmic-ray 
sources. 



4.2. The Correlation Between UHECRs and VCV AGNs 

All sky distribution of AGNs in our selected sample and 
the deflection-corrected positions of the 135 UHECR events 
are presented in Figure |5] We perform the correlation analy- 
sis separately for Auger, AGASA, and Yakutsk events. The 
correlation results for the UHECRs with and without deflec- 
tion corrections are presented in Table [T] We also present 
the correlation results for deflection corrections using the PS 
model ( Prguza & Smida 2003 ) and various GMF components 
and their combinations in our toy model. In Figures |6] and |7] 
we show the results of pair counting and probability analysis 
for the deflection-corrected UHECRs by our GMF toy model. 

In Figures |6] |2l and Table 1, we found that the num- 
ber of Auger UHECR-AGN pairs is more than that of 
simulated isotropic random UHECR samples, which sug- 
gests that the Auger UHECRs are anisotropic and some- 
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<5 n 

Fig. 6. — Number of UHECR-source pairs (dot) for nearby VCV AGNs 
(left column) and Fermi sources (right column), as a function of the angular 
separation (bin size) S. The solid line is the Monte-Carlo simulated average 
number of UHECR-source pairs with errorbar of ilcr, derived from random 
isotropic distribution of simulated cosmic-ray events. 



ho w correlated with VCV AGNs , supporting the results 
of iThe Pierre Auger Collabora tion' ('2007', '2008'). For the 
deflection-corrected Auger UHECRs by our GMF model and 
the PS model, the chance probability is 2 x 10^^, and the 
correlation significance is about 4(7, which is similar to but 
slightly less significant than the results of a chance probabil- 
ity of 1 X 10^^ and correlation significance of 4.7a with- 
out deflection correction. In Table 1, we found that the 
similar correlation results are also presented in the cases 
of various GMF components and the PS model, which in- 
dicate that some of UHECRs probably come from a few 
of AGNs. The marginal correlation significance of a few 
sigma (< 5cr) only suggests that most UHECRs do not 
have AGN counterparts. To further evaluate the deflec- 




6 n 

Fig. 7. — Chance probability P(i5) as a function of the angular separation 5 
of Auger UHECR events from the VCV AGNs and the Fermi sources (dashed 
line), for example. 



tion effect, we use the deflection-corrected Auger events to 
se arch for match pairs, and comp a re the m with the results 
of iThe Pierre Auger Collaboration! (12008) who did not con- 
sider the deflections. Using the 442 nearby VCV AGNs 
from 12th edition, we foun d 15 UHECR- AGN pairs, in- 
stead o f 20 pairs reported by .The Pierre Auger Collaboration! 
(12008!) with the pai-ameters z^ax = 0.017, ip = 3°2, 
and i?th=57 EeV Note that the fewer pairs imply the 
weaker correlation between UHECRs and AGNs after de- 
flection correction. Some of the matc hed pairs reported by 
iThe Pierre Auger Collaboration! (l2008h are probably due to 
random match. The Auger UHECR-AGN correlation is also 
weakened when the update Auger UHECR data are used 
(The Pierre Auger Collaboration 2009). 

As shown in Figures |6]|7] and Table [T] we found no corre- 
lation between the AGASA and Yakutsk events and AGNs. 

4.3. Correlation Between UHECRs and Fermi j-ray Sources 

High energy 7-ray emissions are thought t o be a distinctive 
feature of the possib l e sour ce of UHECRs (iGorbunov et alJ 
l2002h . [Torres et"an (1200 3*) searched for possible correla- 
tion between the fliird EGRET sources and AGASA UHECR 
events, no correlation was found. The GeV AGNs detected 
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by the Fermi/LAT should be most energetic AGNs which 
may have powerful ability to accelerate the charged parti- 
cles to ultrahigh energy bands. Most important is that the 
selected GeV AGNs in the Fermi catalog are relatively com- 
plete above the 7-ray flux threshold. We therefore investi- 
gate the possible correlations between Fermi sources and 
the Auger/AGASAA'akutsk UHECR events considering the 
deflection o f the GMFs. A t the fi nal stage of this work, we 
noticed that lMirabal & Oval ( 1201 Oh did the similar work as us 
and studied the possible correlation between the Fermi/LAT 
First Source Catalog with the public Auger UHECRs. They 
found that 12 of the 27 Auger UHECRs arrived within 3°1 of 
Fermi sources, similar to the matches of artificially random 
distribution of UHECR samples. They concluded for no cor- 
relation. However, the possible correlation of UHECRs and 
Fermi sources is needed to be re-examined after the deflec- 
tions by the GMFs and the GZK effect are considered. 

The ferm//LAT IFGL catalog contains 1451 point sources, 
including 820 identified objects: 62 pulsars, 295 BL Lacs, 278 
FSRQs, 120 other types of AGNs, and 65 other 7-ray sources. 
We remove the Fermi sources with redshift z > 0.024 by 
considering the GZK cutoff, and also neglect the sources lo- 
cated in the Galactic disk with |6| ^ 10°. Figure [8] presents 
the sky distributions of the remaining 635 Fermi/LAT 7-ray 
sources and the deflection-corrected UHECR events. 

The correlation results for the Fermi 7-ray sources and 
UHECRs are listed in Table [T] In the right panel of Figure |6l 
we also present the number of UHECR-source pairs as a func- 
tion of the bin size S. In the case of deflection-corrected Auger 
events by our GMF model, a small excess appears around 
S ^ 0°9, with a chance probabiHty of F ^ 2.5 x lO^^^, and 
a significance level of ^ A.la. A correlation with a chance 
probability of F ^ 6.0 x 10^*, and a significance level of 
^ 3.7(7 is also found for the case with deflection correction 
by the PS model. These correlations with marginal signifi- 
cance probably indicate that a few Fermi sources are related 
to the UHECRs. However, most of the Fermi sources are 
not potential sources of UHECRs. No evidence of correla- 
tion is found in the case of the observed positions of Auger 
events. In the cases of the AGASA and Yakutsk data, no sig- 



nificant correlations are found, though many UHECRs have a 
Fermi source within 3° from their arrival positions. In order 
to evaluate the effect of deflections on the c orrelation anal- 
ysis a nd compare our results with those of iMirabal & Oval 
(I2OIOI) . we also try to use all 1451 sources in the IFGL for 
correlation with Auger UHECRs. We found that the number 
of matched UHECR-source pairs can be reproduced by the 
simulated isotropic UHECR samples for the UHECR sam- 
ple without deflection correcti on, which is consistent w ith the 
conclusion of no corr elation bv lMirabal & Oval (l2010l) . When 
we exclude the sources with |6| ^ 10° and consider the de- 
flection correction by our GMF model, a marginal correlation 
is found with a chance probability of F ^ 9.8 x 10^"^ and a 
significance level of ~ S.lfi. 

Fermi sources contain several types of objects, such as 
pulsars, AGNs, and unidentified sources. It is interesting 
to see if possible correlations exist between some types of 
Fermi sources and UHECRs. The correlation analysis re- 
sults for the Fermi AGNs and UHECRs are presented in Ta- 
ble[T] The Fermi AGNs are weakly correlated with AGASA 
UHECRs after the deflection corrected by our GMF model 
and the PS model. Three hundred seventy-three unidentified 
Fermi sources are also weakly correlated with Auger UHE- 
CRs after the deflection corrected by the PS model. The cor- 
relation results for UHECRs with deflection corrections using 
various GMF components and their combinations in our toy 
model are also presented in Table [T] 

Fermi IFGL catalog has 8 objects of redshift z ^ 0.024: 
NGC 253 (z = 0.001), NGC 4945 (z ^ 0.002), Centaurus 
A (Cen A, z = 0.002), M87 (z = 0.004), ESO 323-G77 
(z = 0.015), NGC 6951 (z = 0.005), NGC 1275 (z = 0.018), 
and M-82 (z = 0.001). We found that 4 of these 8 objects 
have UHECR counterpart(s) within 3° 1 from the m, after the 
deflect ions are corrected by our toy model. Mir abal & Oval 
(120 lOl) found 2 Fermi objects (NGC 4945 and Cen A) within 
3°1 from 3 of the 27 Auger UHECRs. When we consider the 
deflection correction, we find one Auger UHECR event within 
3.1° from NGC 4945, Cen A and ESO 323-G77. NGC 6951 
also has one Yakutsk event within 3°1. 

The Cen A is the nearest FR II radio galaxy dlsraell 



9 



Il998l) . which has bee n long propos e d as a p ossibl e 
source of UHECRs (tCavallol [19781 iRomerol [T996h . 
Cen A was detected at MeV to GeV energ ies b y 
the Fermi/LAT ( The Fermi/LAT Collaboration] |2009|) . 
iThe Pierre Auger Collaboratiom (i2007t i2008i) pointed that 
4 of the 27 events were possibly associated with Cen A 
(e.g. 'Mosk alenko et al.ll2009t iKachelrieB et all 12009). From 
Figure |3j we found that the arrival directions of cosmic ray 
in the region of Cen A are not significantly corrected by the 
GMF model. Two of UHECR events in the 27 pubUshed 
Auger data set are very c lose to Cen A . Considering heavier 
composition of UHECRs. iPiranI (|20 1 Oh suggested that Cen A 
is the only active potential source of heavy nuclei UHECRs 
within a few Mpc for the GZK cutoff. The heavy nuclei 
suffer a larger deflection which can erase any correlation 
with their source. If all detected UHECRs are produced by 
merely Cen A, which come to our Galaxy and suffer different 
deflection via different paths and finally arrive at the Earth 
from various directions, they should show some kind of 
concentration around the source direction for many lighter 
nuclei, depending on the detailed composition and magnetic 
deflection. In Figure |8] there is some indication for such a 
concentration within about 20° near Cen A. While, other 
UHECRs coming from other very different directions may 
have other accelerating sources rather than Cen A. 

NGC 4945 is i dentifi ed as a Seyfert galaxy 
(lY eron-Cettv & Veron' 2006), also known as a s tarburst 
galaxy (Lenc & Tingay 2009). i Lenc & Tingavl (|2009|) identi- 
fied a non-thermal source with a jet-like morphology near the 
AGN of NGC 4945. ESQ 323-G77 i s identified a s a Seyfert 
galaxy having strong Fe/ J emission (lFairalllll986h . however, 
no compact radio core or a radio excess has been detected 
([Corbett et al. 2003). NGC 6951 is known as a LINER galaxy 
and a bipolar outflow which seem s to be associated with 
a nucl ear jet has been reported by IStorchi-Bergmann et al.l 
(l2007l) . However, it is not clear whether the possible radio 
jets from NGC 4945 and NGC 6951 could be the accelerator 
of UHECRs. 

5. SUMMARY AND DISCUSSION 

We collected 135 published UHECR events including 57 
UHECRs recorded by AGASA with energy E > 40 EeV, 
51 events observed by Yakutsk, both located in the north- 
ern hemisphere, and 27 events with energy E ^ 57 EeV 
detected by Auger located in the southern hemisphere. We 
use a new GMF toy model constrained by updated mea- 
surements to evaluate the deflection effects on the arrival 
directions of UHECRs. Considering the possible deflec- 
tion correction by our toy model and the PS model, as well 
as the different magnetic field components in our model, 
we search for the possible correlations of UHECRs with 



nearb y AGNs extracted from the n ew 13th VCV AGN cata- 
log of iVeron-Cetty & Veronl(l2010l) and the Fermi/LAT First 
Source Catalog of 7-ray sources. We found a correlation be- 
tween the Auger UHECR events and nearby VCV AGNs with 
a chance probability of 2 x 10~^, and a significance level of 
^ Aa. Using the same data as lThe Pierre Auger Collaboration! 
( I2OO8I) . we found fewer UHECR-AGN pairs when deflec- 
tion is considered, which implies the weakened correlation. 
A marginal correlation was found between the Auger events 
and the first year Fermi 7-ray sources with a significance 
level of ^ 4cr if the deflection by the GMF model is con- 
sidered. Some Fermi sources of nearby AGNs, NGC 4945, 
ESQ 323-G77, NGC 6951, and Cen A, may be related to 
UHECRs within 3°1. For AGASA and Yakutsk UHECRs, 
no evidence of significant correlation is found for the nearby 
AGNs or the Fermi sources because the matched pairs can be 
reproduced by the simulated random isotropic UHECR sam- 
ples, though some 7-ray point sources are coincident with the 
UHECR events within 2°. 

The correlations of UHECRs with some astrophysical ob- 
j ects suggest that at least some of the UHE CRs are protons 
jThe Pierre Auger Collabo ration 20071 l2008l) . However, most 
UHECRs seem to come from various directions and do not 
associate with known astrophysical objects, which indicates 
that the majority of UHECRs might suffer larger deflections 
in the trajectory, due to either the unknown extragalactic mag- 
netic fields or the heavy nuclei component of UHECRs (PiraQ 
I2OIOI) . The deflection of heavy UHECRs by the GMF models 
is proportional to the charge of nuclei, which leads to a very 
large deflection angle, for example tens of degrees for iron, 
and then any correlation discussed in this work can be dimin- 
ished (Gureev & Troisky 2010). If the primaries of the UHE- 
CRs are heavy nuclei, instead of proton, the identification of 
UHECR sources would be very difficult. Obviously, the un- 
derstanding of UHECR origin will strongly depend on our 
knowledge about the strength and configuration of the Galac- 
tic and extragalactic magnetic fields, which definitely needs 
more measurements (.Haa.2008,) . 
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